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ABSTRACT: The Standard Model of electroweak interactions has had great success in de- 
scribing the observed data over the last three decades. The precision of experimental measure- 
ments affords tests of the Standard Model at the quantum loop level beyond leading order. 
Despite this great success it is important to continue confronting experimental measurements 
with the Standard Model predictions as any deviation would signal new physics. As a fun- 
damental parameter of the Standard Model, the mass of the H^-boson, Mw, is of particular 
importance. Aside from being an important test of the SM itself, a precision measurement of 
Mw can be used to constrain the mass of the Higgs boson, Mh- In this article we review the 
principal experimental techniques for determining Mw and discuss their combination into a 
single precision Mw measurement, which is then used to yield constraints on Mh- We conclude 
by briefly discussing future prospects for precision measurements of the VK-boson mass. 
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1 INTRODUCTION 



The Standard Model of electroweak interactions (SM) theoretically unites the 
electromagnetic and weak forces of nature. It postulates that these forces are 
communicated between the constituent particles of nature, quarks and leptons, 
by carriers known as gauge bosons. In particular, the electromagnetic force is 
carried by the photon, 7, while the weak force is mediated by the neutral Z-boson, 
Z*^, and the charged VT-bosons, W^. As such, the VF-boson is fundamental to 
the Standard Model. Moreover, the mass of the VF-boson, Mw, is a parameter 
of the theory itself, so that a comparison between the experimentally determined 
Mw and the SM prediction provides an important and fundamental test of the 
theory. Alternatively, a precision measurement of M\y can be used to estimate, 
within the framework of the SM, other parameters such as the mass of the Higgs 
boson, Mff- 



1.1 Historical Overview 

The weak force was first inferred from observations of nuclear /3-decay, n 
p+e~ +Ve—. In 1935, Fermi postulated the first theory of weak interactions. The 
form of the interaction was taken to be analogous with that of the electromagnetic 
interaction, and was characterized by a "coupling" (or strength) parameter — the 
Fermi constant, Gp- By comparing interaction rates, the strength of the weak 
force was estimated to be about 10~^ that of the electromagnetic force. Fermi's 
theory very successfully described low energy weak interactions, but violated 
unitarity at high energy. 

In 1967 Glashow, Weinberg and Salam proposed the electroweak SU{2) x U{1) 
gauge theory, which unifies the weak and electromagnetic forces Q. The theory 
postulated that the weak force was mediated by massive particles, the W- and 
Z-bosons, and predicted their masses to be of order 10^ GeV Q The discovery 
of the VF-boson in 1983, with a mass of 81 it 5 GeV Q, was a great success for 
the electroweak theory. More rigorous tests of the theory required more precise 
determinations of the boson masses. 

Over the past 15 years a variety of experiments have measured the mass of 
the VF-boson with ever improving precision. The first measurements were made 
at the CERN SppS collider @ by the UAl |] and UA2 |] experiments. The 
UA2 experiment made the first measurement of the PF-boson mass at a relative 
precision below 1% |p. The CDF and D0 experiments at the Fermilab 
TeVatron Q, another pp collider, were the first to push the precision to the 
0.1% level. More recently, measurements made at the CERN e+e~ -collider, LEP, 
by the Aleph |1^, Delphi [|11|, L3 Q, and Opal Q experiments, have also 



reached relative precisions of 0.1%. Combining all these measurements yields a 
relative precision of 0.05% and affords stringent tests of the Standard Model. 
In particular, due to radiative corrections, such precision measurements offer 
indirect constraints on the mass of the Higgs boson. 



1.2 The Electroweak Theory 

In the SU{2) x ^7(1) electroweak theory, local gauge invariance is achieved by 
introducing four massless bosons, an isovector triplet = (Wg^, VFf , 1^2^), 
and an isosinglet, Bq. In analogy to the electromagnetic case, the electroweak 

^Here and throughout this article we use units of ?i = c = 1 
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Lagrangian can be expressed as a product of currents and coupling parameters: 

C = gJ^-W^^ + g'4B^ (1) 

where and Jy are the weak isospin and hypercharge currents of the physical 
fermions {i.e. quarks and leptons), respectively, and g and g' are their couplings 
to the and Bq fields. The weak quantum numbers are related to the electric 
charge, Q,hy Q = I3 + Y/2, where I3 is the third component of the weak isospin 
associated with the SU{2) group and Y is the weak hypercharge associated with 
the U{1) group ||l^. The fact that the associated bosons are massless implies that 
the weak field is a long-range (infinite) field, in contradiction with experimental 
evidence. This short-coming can be addressed by imparting mass to the vector 
bosons, which is achieved by spontaneously breaking the SU{2) x L'^(l) symmetry 
with the introduction of an additional field. Demanding that the theory be valid 
to high energies and remain renormalizable, a necessary condition in order to 
extract meaningful theoretical predictions, constrains the form of this additional 
field. The simplest solution introduces a complex scalar isodoublet, the Higgs 



field, with one component having a vacuum expectation value v > [15|. The 
physical boson fields can then be expressed as 

f cos 6w sin 9\Y \ f Bq 



for the charged VF-bosons, W^, the neutral Z-boson, , and the photon, AP , 
respectively. The weak mixing angle, 0w, relates the SU{2) and U{1) coupling 
constants to the electromagnetic coupling constant {i.e. the fine structure con- 
stant), a, by 

g^ = Ana/ sin^ Oy/-, 5'^ = Aixa/ cos^ Oy/- (3) 
The gauge boson masses are given by 



Mw = gv/2, Mz = v^g^ + g'^/2, Ma = (4) 

corresponding to the massive Wi: and Z^-bosons and the massless photon, 
respectively. Equations ^ and Q yield the following relationship, 
sin2 9w = l- {Mw/Mzf. 

At low energies, the electroweak theory is equivalent to the Fermi theory of 
weak interactions. Comparing the electroweak Lagrangian in Equation |^ to 
Fermi's expression for the weak interaction yields the following equality, 
Gf = g^l{A^Mw^) = 7rQ/(^/2Mvy2sin2^vy)- This can be rewritten as 

relating the mass of the M/^-boson, the mass of the Z-boson, the fine structure 
constant and the Fermi constant, so that a measure of three yields a prediction of 
the fourth. To obtain theoretical predictions of an precision comparable to that of 
the experimental determinations of these parameters, radiative corrections must 
be included. These corrections can be incorporated by rewriting Equation |^ as: 

2 { ^ Mw'^ \ vra f 1 
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where the effects of the radiative corrections are included in the additional term, 
Ar. The corrections can be separated into three main pieces, 

Ar = Aa + Ap{mto/) + Ax{HMh/Mz)), (7) 

which include the running of the fine structure constant, Aa, a quadratic depen- 
dence on the top quark mass, Ap, and a logarithmic dependence on the mass of 



the Higgs boson, A^ |16|. This last dependence is a unique consequence of the 
non-Abelian gauge structure of the electroweak theory, which allows interactions 
among the gauge bosons themselves. It is because of these radiative corrections 
that precision measurements of Gp, ct, Mz and M\y when compared with theo- 
retical calculations, can yield constraints on mtop and Mh |17]. 



2 MEASUREMENTS OF Mw AT pp COLLIDERS 

2.1 Measurement Techniques 
2.1.1 1^-BOSON PRODUCTION 

There have been two pp colliders with sufficient center-of-mass energy (y/s) to 
produce Il^-bosons, the SppS at CERN (y^=630 GeV) and the Tevatron at 
Fermilab (y^=1.8 TeV). Figure]^ shows the most important subprocesses for W- 
boson production in pp collisions. At these center-of-mass energies, the dominant 
subprocess is qq W. Z-bosons, which form an essential control sample, are 
produced via analogous processes. 




Figure 1: Feynman Diagrams for Il^-boson production in pp collisions. 



The VF-boson mass measurements from these colliders all make use of the 
W ^ eV and W — > /uz7 decay channels. Electrons and muonsQ are easy to trigger 
on and their momenta can be measured very precisely. Moreover, W and Z- 
bosons are the dominant source of isolated, high-py electrons and muons in pp 
collisions. Therefore, samples of W and Z-decays involving electrons and muons 
can be identified with very little background. Purely hadronic decays of the W- 
boson are swamped by QCD background. Decays involving r leptons are difficult 
to identify because the r leptons decay before they enter the detector. 

The cross sections for W and Z production in pp collisions are large, axB = 680 
pb at ^=630 GeV ^ and 2.3 nb at y/s=1.8 TeV |l| for M^-bosons, where B 
is the leptonic branching fraction. For Z-bosons the corresponding values are 
about 10 times smaller. 

In the following, we refer to a coordinate system which has its origin at the 
average pp-collision point. The z-axis is defined by the proton beam. The y-axis 

^Here and throughout this article charge conjugation is imphed. 
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points up. The x-axis points along the horizontal. Since the parton center-of-mass 
frame is boosted along the beam direction, momentum components transverse to 
the beam are especially important. They are are denoted by a subscript t- The 
beams are unpolarized so that there is an inherent azimuthal symmetry. Thus 
it is often convenient to work in a cylindrical coordinate system in which (f) is 
the angle with the x-axis in the x-y plane. The longitudinal phase space is most 
conveniently expressed in terms of the pseudo-rapidity, rj = — In tan {9/2), which 
is related to the polar angle 9. 

The detectors have approximate azimuthal and forward-backward symmetry. 
They are constructed to cover as large a region in pseudo-rapidity as possible. 
Prom inside out, they typically consist of several sub-detectors: a tracking system 
to measure the trajectories of charged particles; a calorimeter to measure the 
energy of electrons, photons, and hadrons; and a muon detection system. The 
tracking system may be located in a magnetic field to determine the momentum 
of charged particles from the curvature of their trajectories. 

2.1.2 EVENT CHARACTERISTICS 

The detectors register the charged lepton from the decay of the M^-boson, while 
the neutrino escapes without detection. The initial proton and antiproton break 
up in the collision and the fragments hadronizc, contributing additional particles 
to the event. The hadronization of final-state quarks or gluons also contributes 
particles, which may form jets if the initial parton had sufficient transverse mo- 
mentum. We refer to all particles, except the VF-boson and its decay products, 
as the underlying event. Some of the particles of the underlying event escape 
through the beam pipe and are not detected at all. These particles may carry 
a substantial momentum component along the beam axis, but they carry little 
momentum transverse to the beam. 

The transverse momenta of all final-state particles must add to zero, because 
the initial pp momentum is zero and momentum is conserved. Since the unde- 
tected neutrino carries away substantial momentum, the transverse momenta of 
all observed final-state particles do not add to zero. The apparent transverse 
momentum imbalance is called "missing py" • 

The particles of the underlying event that fall within the detector acceptance 
cannot all be detected individually. The detector measures the sum of the energy 
of all particles incident on one calorimeter segment. The quantity 



gives an approximate measurement of the total transverse momentum of all un- 
derlying event particles. The sum runs over all calorimeter cells, except those 
assigned to the charged lepton. is the energy in cell i. The unit vector i forms 
a right angle with the beam axis and points from the beam to cell i. 

Thus, the basic obscrvables are the momentum of the charged lepton (j/, ^=e or 
fi) and the sum of the transverse momenta of the particles in the underlying event 
{ut) which we call the recoil momentum. From these, the transverse momenta of 
the VF-boson {p^ = —ut) and the neutrino = —ut — p^) can be inferred. A 
high-pT charged lepton and large missing pT form the characteristic signature of 
H^-boson decay events. Z-decay events are characterized by two charged leptons 
with high pt- There are no high-pT neutrinos in such Z-decays and therefore no 
significant missing px is expected. 




(8) 
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2.1.3 MASS MEASUREMENT METHOD 

It is not possible to reconstruct the invariant mass of the VF-boson because there 
is no measurement of the momentum component of the neutrino along the beam 
axis. In addition, the VF-bosons are neither produced at rest nor arc they the 
only particles produced in the collisions. Therefore, a precision measurement of 
the I^-boson mass using PF-bosons produced in pp collisions poses a particular 
challenge. 

The most precise measurements of the mass of the W^-boson are based on the 
transverse mass of the charged lepton-neutrino pair, defined as 

rriT = Y^2pfp^ (1 - cos {(j/ - (j)")). (9) 

The advantage of the transverse mass is its invariance under boosts along the 
beam axis. Boosts transverse to the beam axis only give rise to corrections 
of order {p!^ /E^)'^. On the other hand, the transverse mass depends on the 
measurement of the recoil momentum ut and all the associated systematic effects. 

An alternative method to determine My\r uses the pT spectrum of the lepton. 
This has the advantage of being insensitive to ut- However, it is affected by 
the boost of the W-boson transverse to the beam axis to order /E^ and is 
therefore much more sensitive to systematic effects associated with the production 
of the ly-bosons. 

In principle, the charged lepton momentum or the transverse momentum of the 
neutrino can also be used to measure the W-^-boson mass. However the former is 
sensitive to boosts in all directions and the latter suffers from poor resolution. 
These variables serve mainly as cross checks. 

It is not possible to describe the spectra of the variables mentioned above 
analytically. They have to be calculated numerically using a Monte Carlo model 
that takes into account the mechanisms for production and decay of VF-bosons, 
and detector effects. 

2.1.4 BACKGROUNDS 

The backgrounds to the W ^ £u signal are W ^ tV —>■ iUUu (1-2%), hadronic 
backgrounds (1% for W eV, <1% for W fiV), Z ^ £+£' (<1% for W eV, 
4% for W i-w), and cosmic rays (<Cl%). Hadronic backgrounds arise hadrons, 
that fake the charged lepton signature. Z f.'^fr decays can enter the W 
sample if one of the Icptons escapes detection. The quoted percentages give the 
approximate residual fractions of background events in the final W samples. The 
precise background contamination depends on the details of the event selection 
and the detector. They have to be taken into account in the measurement to avoid 
biasing the result. The normalization and shapes of the background spectra are 
determined from control data samples and detailed Monte Carlo simulations. 

2.1.5 EVENT SELECTION 

The event selection consists of the identification of the charged lepton and a set 
of kinematic and topological cuts. The selection criteria have to achieve two com- 
peting goals: to reject backgrounds efficiently and to minimize any biases to the 
selected event sample. Kinematic cuts, requiring the momentum of the charged 
lepton and missing pT above a threshold (typically 25 GeV), are easy to simu- 
late and reduce backgrounds significantly. VF-bosons with very high transverse 
momenta do not add to the statistical significance of the mass measurement. 
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because their transverse mass and lepton pT spectra are very broad and carry 
little mass information. In addition, their recoil response is difficult to simulate 
and they are subject to higher background contamination. Thus, such events are 
usually eliminated from the sample by requiring that the P^-boson is below 
some threshold and/or that there are no high-py jets in the events. 

An electron is typically identified as an energy deposit in the calorimeter, con- 
sistent in shape with an electromagnetic shower, and a track that points to it. 
Since these electrons are highly relativistic, their momenta can be equated to the 
energy measured in the calorimeter. A muon is typically identified as a track stub 
in the muon detection system that matches a track in the tracking system and 
energy deposits in the calorimeter, small enough to be consistent with the pas- 
sage of a minimum-ionizing particle. These criteria reduce contamination from 
hadronic backgrounds. However, both criteria inherently require the lepton to be 
isolated from other activity in the event. This biases the selection towards event 
topologies in which the charged lepton is emitted along the direction of motion of 
the W-boson. In such events, the probability is smaller that the lepton overlaps 
with a recoil particle. Since the boost of the VF-boson increases the lepton pT 
in the lab frame, these events tend to have harder charged leptons and softer 
neutrinos. This bias docs not affect the transverse mass spectrum significantly, 
but it must be understood to predict the pT spectra of the charged leptons and 
the neutrinos correctly. 

Specific cuts are required to reject events due to an accidental coincidence 
between a cosmic ray traversing the detector and a beam crossing. 

2.1.6 MONTE CARLO MODEL 

In this section a generic description of the Monte Carlo models is given. The 

sections below, describing the individual measurements, highlight significant ex- 
periment specific deviations. To keep statistical fluctuations in the Monte Carlo 
simulation negligible, many millions of 1^-decay events have to be generated. To 
simulate such large event samples, parameterized Monte Carlo algorithms for TV- 
boson production and decay, and detector modeling were developed specifically 
for the VF-mass measurements. 

First, the VF-bosons arc generated. Their px distribution is determined theo- 
retically from QCD-based calculations, empirically from the observed px distri- 
bution of .Z^-bosons, or by a combination of both. The rapidity distribution of 
the generated IF-bosons depends on the momentum distribution of the partons 
inside the proton. To determine it, a specific set of parton distribution functions 
must be chosen. The mass distribution of the generated PF-bosons is a relativistic 
Breit-Wigner curve with peak at the hypothesized value of the VF-boson mass and 
s-dependent width, given by the Standard Model expectation. This mass spec- 
trum is skewed towards lower mass values due to the momentum distribution of 
the partons inside the proton. 

Next, the VF-bosons are allowed to decay. At lowest order, the angular dis- 
tribution of the charged leptons is da/dcosO* oc (1 — ^qcosO*)'^, where 6* is the 
scattering angle of the charged lepton in the rest frame of the VF-boson, q the 
charge of the lepton, and ^ the helicity of the 1/F-boson. In most events the initial 
quark comes from the proton and ^ equals —1. In the much less likely case that 
the initial antiquark comes from the proton ^ equals -|-1. Higher-order QCD pro- 
cesses modify the angular distribution of the leptons. Radiative decays, in which 
W ^7, modify the momentum spectrum of the leptons. The Monte Carlo 
models either include these effects or corrections are applied to the results. 
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The decay W ^ tV ^ IWv leads to events, that are topologically indistin- 
guishable from W tu. These can be calculated precisely in the framework of 
the Standard Model and are typically included in the Monte Carlo model. 

Finally, the Monte Carlo model must account for detector effects. The simula- 
tion has now generated the "true" momenta of the H^-boson (and thus the recoil 
momentum) and the charged lepton in the event. These are modified to account 
for experimental resolutions, biases, and efficiencies. Adding random Gaussian 
uncertainties to the observables simulates resolution effects. The widths of these 
Gaussian distributions are parameterized in detector-specific ways. Other effects 
accounted for include the response of the detector to the charged lepton and 
to the underlying event. Also modeled are the imperfect separation of energy 
deposits between lepton and underlying event, which leads to biases in lepton 
and recoil momentum measurements, and selection efficiencies that depend on 
kinematics or topology of the events. 

Events due to the process pp ^ Z + X, Z ^ i^i^ constitute an extremely 
important control sample to determine these effects. Comparing the observed Z 
peak to the known Z-boson mass calibrates the energy or momentum response 
to charged leptons. The observed width of the Z peak measures the energy or 
momentum resolution for charged leptons. The Z-boson px can be measured 
directly using the charged leptons from its decay and indirectly from the recoil 
momentum. By comparing both determinations, the Z events also serve to cali- 
brate the recoil momentum response of the detector relative to the charged lepton 
response. 

2.1.7 MASS MEASUREMENT 

The Monte Carlo model predicts the shape of the transverse mass and the lepton 
Pt spectra from W ^ tu decays as a function of the hypothesized value of 
the VF-boson mass. These are added to the estimated background spectra and 
normalized to obtain probability density functions for a maximum likelihood fit 
to the spectra from the collider data. The statistical uncertainty in the fit to the 
niT spectrum is typically 11 GeV/(number of events)^/^. 

Figures |2| and ^ show representative spectra of transverse mass and lepton pT- 
These spectra were measured by the D0 experiment. The points indicate the 
collider data, the line indicates the Monte Carlo prediction that agrees best with 
the data, and the shaded region indicates the estimated background contribution. 

2.1.8 SYSTEMATIC UNCERTAINTIES 

All inputs of the Monte Carlo model contribute to the systematic uncertainty 
of the measurement. Each such contribution is estimated by varying the input 
within its 68% confidence interval to determine the resulting change in the mea- 
sured VF-boson mass. The total systematic uncertainty is the sum in quadrature 
of all such contributions. Most of the model parameters are constrained by con- 
trol data samples, most notably by the Z — > samples. In most cases the 
precision with which these parameters can be determined is limited by the size 
of the control samples, so that these uncertainties are really statistical in nature. 
This means that they can be quantified in a well-defined way. There are some 
cases in which no rigorous confidence interval can be defined, as is usually more 
typical of systematic uncertainties. 

The following paragraphs elaborate on the most important categories of sys- 
tematic uncertainties. The values quoted for each uncertainty are typical of the 
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Figure 2: Transverse-mass spectrum from VF-decays measured by D0 [p^. 




10 



Glenzinski & Heintz 



measurements from the Tevatron using a fit to the niT spectrum from a data 
sample of about 100 pb~^. 

Lepton energy /momentum scale: (70-85 MeV) This is the most important 
systematic effect. Since ah detector responses are cahbrated against the 
charged leptons using the Z sample, the measured W mass simply scales 
with the lepton scale. It can be set in two ways. 

One method is to calibrate the lepton scale so that the Z e^e" and Z 
jJi^ lJi~ mass peaks (Fig. ^ agree with the known Z-boson mass ||2l|. This has 
the advantage that the uncertainty is dominated by statistical fluctuations 
in the Z sample, «3 GeV/(number of events)^/^, and little extrapolation 
is needed to the energies of leptons from VF-decays. Uncertainties in the 
extrapolation can be limited by using other resonances, such as J /ip e~^e~ 
or TT^ ^77- If the scale calibration of the charged leptons is tied to the Z- 
boson mass, the measured quantity is really the ratio of the W and Z-boson 
masses, rather than the VF-boson mass. Given that the Z-boson mass is 
known so precisely, the two quantities are of course de facto equivalent. 




Detectors with a magnetic tracking system can alternatively calibrate the 
momentum measurement for charged tracks, using for example J/V' — 
decays (Fig. ||), and then extrapolate to the momentum of leptons 
from M^-decays. The calorimeter must then be calibrated against the track 
momentum using the ratios of energy and momentum (E/p) measured for 
electrons from I^-decays. The advantage of the latter method lies in the 
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higher precision of the track momentum cahbration. Its disadvantage are 
the systematic effects associated with the extrapolation to higher momenta 
and the effects of radiation on the E/p spectrum. 



J/if/ — > /x/j. Data 
Simulation 




3000 3020 3040 3060 3080 3100 3120 3140 3160 3180 3200 



m 



mi 



(MeV) 



Figure 5: Mass spectrum from J ftp n'^fj, decays measured by CDF |22]. 



Lepton energy /momentum resolution: (20-25 MeV) The electron energy res- 
olution can be modeled as a/E = {S'^/E -|-C^)^/^, where E is the electron 
energy, S the sampling term and C the constant term. The value of S is 
taken from beam tests and C is chosen so that the width of the Z peak pre- 
dicted by the Monte Carlo model agrees with the Z peak from collider data. 
For muons, the transverse momentum resolution is of the form a/pj^ = k, 
where pT is the transverse momentum of the muon and k a constant chosen 
to match the widths of the Z peaks from the Monte Carlo model and the 
collider data. 

Recoil model: (30-40 MeV) This uncertainty arises from the parameters de- 
scribing the response and resolution of the detector to the underlying event. 
These are determined from Z-decays and to a lesser extent from VK-decays. 

Lepton removal: (~ 15 MeV) This describes uncertainties in the corrections 
to the recoil momentum for the imperfect separation of energy deposits 
between the charged lepton and underlying event. Some particles from the 
underlying event inevitably overlap with the charged lepton in the calorime- 
ter. Their energies are not included in the calculation of ut- The correction 
is equal to the average energy deposited by the underlying event in an ap- 
propriately sized calorimeter segment in the W data sample. 

Proton structure: (10-20 MeV) This uncertainty characterizes the variations in 
the result between different choices of parton distribution functions. While 
relatively small, this uncertainty is completely correlated for all measure- 
ments at pp colliders. For any given set of parton distribution functions, the 
variation in the measured VF-boson mass is strongly correlated with the vari- 
ation in the predicted forward-backward charge asymmetry m.pp ^ tu 
[p2| . Recently, increasingly precise measurements of this asymmetry by the 
CDF collaboration P3|, pj] have helped constrain the parton distribution 



functions (e.g. |25, |26|) and reduce the resulting uncertainty in the W^-boson 
mass measurement. The agreement of recent parton distribution functions 
with the measured asymmetry is shown in Fig. n. Since no complete er- 
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ror matrices are available for parton distribution functions, this uncertainty 
cannot be evaluated in a statistically rigorous fashion. 
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W pt spectrum: (~ 10 MeV) The shape of the W pT distribution can be con- 
strained by theoretical calculations in conjunction with the observed Z pT 
distribution. For low px, the W px spectrum cannot be calculated perturba- 
tively. One approach [27, ^ uses the Collins-Soper-Sterman resummation 
formalism |3^ which contains universal empirical parameters. These 
parameters can be constrained by requiring the calculation to agree with 
the observed Z pT spectrum. Another, essentially equivalent, approach is 
to use the observed Z px spectrum directly and convert it to a prediction 
for the W pT spectrum using the ratio of the calculated W and Z pT spec- 
tra. The precision of both approaches is limited by statistical fluctuations 
in the Z data sample. The measured W px spectrum does not provide a 
very stringent constraint because its shape at low pT is dominated by the 
recoil response of the detector. The Z px spectrum, on the other hand, 
can be measured independently of the recoil response using the two charged 
leptons from the decay of the Z. 

Higher-order corrections: (10-20 MeV) Calculations of radiative H^-boson de- 
cays involving one photon |3^ and two photons |^] are available. The 
precision of these corrections is limited by experimental factors which de- 
termine whether the photons and the lepton are resolved. 

Backgrounds: (5-25 MeV) The hadronic background dominates this uncer- 
tainty for the — > ez7 channel. Normalization and shape are determined 
from control data samples. The Z — > background dominates the 

W fiV channel. Uncertainties in the tracking efficiency at high \ri\ and in 
the parton distribution functions give rise to this uncertainty. 
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2.2 Individual Measurements 

2.2.1 UA2 EXPERIMENT 

The UA2 collaboration published the first measurement of the W^-boson mass 
with a precision below 1% |^]. This was superceded by an improved result |^] 
based on 13 pb~^ of data taken in 1988-1990 at the CERN SppS collider at 
^=630 GeV. 

The UA2 detector consists of a calorimeter which covers the pseudo-rapidity 
range |r/| < 3. It consists of lead and iron absorber plates interspersed with 
scintillators and wavelength shifter readout. The electromagnetic section is 17- 
24 radiation lengths deep and segmented into elements covering 15° in azimuth 
and approximately 0.2 units in pseudo-rapidity. The electron energy resolution 
is a/E = 17%/ y^E/GeV. The hadronic section is four interaction lengths deep. 
Inside the calorimeter are nested cylindrical tracking detectors. From inside out, 
they are: a drift chamber with arrays of silicon pad counters on either side, a 
transition radiation detector, and a scintillating fiber detector. The detector has 
no magnetic field. 

The VF-boson mass measurement uses the W ^ eV and Z — > e^e~ decay chan- 
nels. The selection for the W-event sample requires an electron in the central 
calorimeter, >20 GeV, >20 GeV, ut < 20 GeV, and AO < mr < 120 GeV. 
The Monte Carlo model calculates the VF-boson rapidity from HRMSB structure 
functions |35|. The spectrum of is taken from a calculation |28|, modified by 
an empirical distortion function. The distortion function is determined by com- 
paring the spectrum of p^ predicted by the same calculation with the observed 
Pj, distribution. The recoil response model has two parameters: resolution (de- 
pendent on the total energy measured in the event), and offset (dependent on 
Px'). Both were tuned using the Z — > e^e~ sample and requiring that the mean 
predicted by the model agrees with the data. 

Two Z samples are used. Sample 1 requires two central electrons, which must 
be inside the fiducial volume of the calorimeter within \r]\ < 0.8. Sample 2 requires 
one central electron and one electron outside the central acceptance region. The 
energy of the "outside" electron is rescaled so that all transverse momentum 
components along the outer bisector of the two electron directions add to zero. 

A fit to the transverse mass spectrum gives 80.84±0.22±0.83 GeV0. The Z mass 
is measured to be 91.74 it 0.28 ±0.93 GeV using both Z samples. In all cases, the 
systematic uncertainties are dominated by the electron energy scale calibration. 
In the ratio Mw /Mz the energy scale and other systematic uncertainties partially 
cancel. UA2 finds Mh^/Mz=0.8813±0.0036±0.0019. Using the current Z mass 
of 91.187±0.002 GeV |l| gives ^ Mw = 80.36 ± 0.33 ± 0.17 GeV. 

2.2.2 CDF EXPERIMENT 

The CDF collaboration has performed measurements of the VF-boson mass using 
data sets from three running periods of the Fermilab Tevatron: 1988/89 



1992/93 HI, and 1994-96 A pubhcation of the results from the 1994-96 

data is in preparation. 

The CDF detector ^ is a multipurpose magnetic spectrometer. Tracking 
detectors are surrounded by a solenoid, that provides an axial magnetic field of 1.4 
T. The vertex time-projection chamber measures the position of the pp-collision 

^Whenever two uncertainties are given, the first is due to statistical fluctuation, the second 
to systematic effects. 
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point along the z-axis with 1 mm resolution. The central tracking chamber has 84 
layers of wires and covers 40° < 9 < 140°. The transverse momentum resolution 
is a/pl = 0.0011/GeV. 

The central calorimeter covers \rj\ < 1.1. The electromagnetic section consists 
of lead plates interleaved with scintillator. Including chamber wall and solenoid, 
it is 19 radiation lengths deep and segmented into projective towers covering 
A(j) X Ar/=15°x0.1. The electron energy resolution is o" = 13.5%y^E sinO/GeV. 
Proportional chambers after 6 radiation lengths measure the shower centroid 
position to 3 mm. The hadron calorimeter is made of iron-scintillator shower 
counters. Outside the central region (1.1 < |r/| < 4.2) the calorimeter is made of 
gas proportional chambers with cathode pad readout. 

Muon chambers are located 3.5 m from the beam behind 5 nuclear absorption 
lengths and cover |r/| < 0.6. 

CDF use both the W ^ eV and W channels. Events are selected with 

Pj. > 25 GeV and pj. > 25 GeV. There must be no high-p-j- tracks or energetic 
clusters in the calorimeter in addition to the charged lepton. For the 1992/93 
data sample, ut < 20 GeV is required. 

The muon momentum scale is based on a calibration of the tracking system 
to the J/tp mass. The electron energy scale is set using E/p for electrons from 
VK-decays. The calibration is checked using the Z mass from Z — > e''~e~ decays, 
91.12±0.52 GeV (1988/89) and 90.88 ± 0.19 ± 0.20 GeV (1992/93), using the 
same calibration as for W ^ eV events. 

For the analysis of the 1994-96 data, the electron energy scale determined 
by the E/p technique results in a Z-mass peak from Z e~^e~ decays 3.9 
standard deviations below the known Z mass. Thus this technique is not used to 
determine the W mass. Instead, the muon momentum and electron energy scales 
are calibrated using the observed Z-mass peaks. 

For the analysis of the 1988/89 data, the Monte Carlo model uses MRS-B par- 
ton distribution functions [^9| as the nominal choice. The transverse momentum 
distribution of the T^-bosons is obtained from the observed p^ distribution by an 
unfolding procedure. The results from the fits to the rriT spectra in both decay 
channels are listed in Table |. Both combined give Mw = 79.91 ± 0.39 GeV. 

For the analysis of the 1992/93 data, parton distribution functions are re- 
stricted to those consistent with the measured charge asymmetry va. W ^ tU 
decays |^3|. The Monte Carlo model uses MRSD— ' [^] as the nominal choice. 
The transverse momentum distribution of the VF-bosons is obtained from the 
observed pj< distribution, corrected for electron resolution and scaled so that the 
spectrum of the component of ut perpendicular to the direction of the charged 
lepton agrees with the W data. The underlying event model uses a lookup table 
of Ut versus generated p^, built from the Z — > e^e~ event sample. The results 
from the fits to the rriT spectra in both decay channels are listed in Table ||. They 
combine to Mw = 80.41 ± 0.18 GeV. 

For the analysis of the 1994/95 data, MRS-R2 |41| parton distribution functions 
are used. The p^ spectrum is derived from the observed p^ spectrum, corrected 



based on a theoretical calculation of the ratio of the p^ and p^ spectra |28|| . 
The parameterized recoil model is tuned to W and Z data. The results from 
the fits to the rriT spectra in both decay channels are listed in Table |l|. Their 
combined value is Mw = 80.470 ± 0.089 GeV. 

AH CDF measurements combined give Mw = 80.433 ± 0.079 GeV. 
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2.2.3 D0 EXPERIMENT 

The D0 collaboration has published three measurements of the T^-boson mass 
using the W ^ eV channel. Two measurements, using data from 1992/93 p^ ] 
and 1994-96 [20|, use electrons in the central calorimeter. The third uses data 



with the electron in the end calorimeters |42]. 

The D0 detector |^] is based on a hermetic uranium-liquid argon sampling 
calorimeter, which encloses a non-magnetic tracking system and is surrounded 
by a muon spectrometer. 

The tracking system consists of nested cylindrical sub-detectors: a vertex drift 
chamber, a transition radiation detector, and a central drift chamber, covering the 
pseudo-rapidity region \r]\ < 1. Forward drift chambers on either side extend the 
tracking coverage to |r/| < 3. The chambers provide measurements of direction 
and energy loss of charged particles. 

The calorimeter is housed in three cryostats. The central calorimeter (CC) 
covers \r]\ < 1 and the two end calorimeters (EC) cover 1 < |r/| < 4. The elec- 
tromagnetic section is 21 radiation lengths deep and segmented radially into four 
layers and laterally into towers covering A(/> x Ary=0. 1x0.1. It measures the en- 
ergy of electromagnetic showers with a resolution of a/E = 13.5% /\/E sin 6 / GeV 
and the shower centroid position with a resolution of 2.5 mm in azimuthal direc- 
tion. The hadron calorimeter is 7-9 nuclear interaction lengths deep and provides 
hermetic coverage without projective cracks. 

The event selection for W-decay events requires pj, > 25 GeV, pj< > 25 GeV, 
and ut < 15 GeV. 

The ly-boson pT and rapidity spectra are determined by a theoretical calcula- 



tion 1 27], constrained against the observed pj^ spectrum, and the MRS A' parton 



distribution functions |25] for the 1992/93 data, and the MRST parton distribu- 
tion functions |44] for the 1994-96 data. The electron energy scale calibration 
is mainly based on the observed Z peak. The energy spread of electrons from 
Z-decays, and signals from J/tp and vr'^ decays limit nonlinearities. Electron 
resolution and the recoil model parameters are determined from the Z data. 

Based on the fit to the niT spectrum from the 1992/93 data, the D0 collab- 
oration measure^ Mw = 80.35 ± 0.21 ± 0.15 GeV. Based on the 1994-96 data, 
the D0 collaboration measures the VF-boson mass using the m-p, pq., and pj^ 
spectra for electrons in CC and EC. Using the complete 6x6 covariance matrix, 
these results are combined to Mw = 80.498 ± 0.095 GeV with = 5.1 for five 
degrees of freedom. By increasing the acceptance for electrons to pseudo-rapidity 
between —2.5 and 2.5, the sensitivity to the rapidity spectrum of the M^-bosons 
is greatly reduced. This is reflected in the reduced uncertainty due to proton 
structure. This uncertainty is 15 MeV if only central electrons are included and 
7 MeV if also electrons in EC are accepted. 

All D0 measurements combined give Mw = 80.482 ± 0.091 GeV. 



2.3 Combination of pp Collider Results 

Table |l| lists the individual measurements for comparison in sequence of their 
publication. The number of VF-boson events given reflect the number of events 
included in the fit to the transverse mass spectrum. The number of Z events 
is given if the Z data were used to calibrate the lepton scale. The statistical 
uncertainty reflects statistical fluctuations in the W data sample. The scale 
uncertainty refers to the uncertainty in the lepton momentum scale calibration. 

^Updated uncertainties [B0[ relative to original publication. 
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If the Z data are used to calibrate the lepton scale, this uncertainty is dominated 
by statistical fluctuations in the Z data sample. The systematic uncertainty 
reflects all other systematic effects. All uncertainties are rounded to the nearest 
5 MeV. 



Table 1: Comparison of individual Mw measurements from pp colliders 



measurement 



events 



M\Y stat scale syst total 





W 


Z 


(GeV) 




(MeV) 




eU 


2065 


251 


80.36 


220 


260 


150 


370 


eV 


1130 


N/A 


79.91 


350 


190 


240 


465 


fiV 


592 


N/A 


79.90 


530 


80 


315 


620 


eV 


5718 


N/A 


80.49 


145 


120 


130 


230 


fjjy 


3268 


N/A 


80.31 


205 


50 


120 


245 


eV 


5982 


366 


80.35 


140 


160 


145 


255 


eV 


23068 


2179 


80.44 


70 


70 


60 


115 


eV 


11090 


1687 


80.69 


110 


190 


75 


230 


eV 


30100 


1600 


80.47 


65 


75 


55 


115 




14700 


1800 


80.47 


100 


85 


55 


145 




CDF||] 

D0 g 

D0 [20] 
D0 [43] 
CDF[37] 



In combining the results from the three pp-collider experiments, correlations 
must be accounted for. Since the Monte Carlo models used by the three experi- 
ments were tuned independently based on experimental data, the detector models 
are certainly independent. The uncertainties due to higher order corrections are 
dominated by independent experimental uncertainties. The constraints on the 
W pt spectra are dominated by statistical fluctuations in the respective Z data 
samples and are therefore uncorrelated as well. Thus, the only significant cor- 
relation originates from the common uncertainty in the structure of the proton. 
Table ^ summarizes the combined data. 



Table 2: Summary of combined measurements of the VF-boson mass at pp colliders 







uncertainty (MeV) 


experiment 


Mw (GeV) 


total correlated 


UA2 


80.363 


371 85 


CDF 


80.433 


79 25 


D0 


80.482 


91 9 



The individual results combine to 

MwiPp) = 80.452 ± 0.060 GeV 

with = 0.23. 



(10) 
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Figure 7: The tree level diagrams for the process e~^e 
neutrino exchange, and s-channel 7 and Z exchange. 




W^W~: t-channel 



3 MEASUREMENTS OF Mw AT LEP 

From 1989-1995 the Large Electron-Positron collider (LEP) at CERN provided 
e~^e~ collisions at center-of-mass energies at or near the Z-boson mass. Since 
1996, LEP has been running at center-of-mass energies above the Vl^-pair pro- 
duction threshold, ^/s > 2Mv/. The collider provides data to four experiments, 
Aleph, Delphi, L3, and Opal. While the LEPl program afforded precision 
measurements of the Z-boson mass, the LEP2 program provides the opportunity 
to precisely measure the VK-boson mass]^. 

At LEP2 energies VK-bosons are predominantly produced in pairs through the 
reaction e"*"e~ — > W^W~ , whose tree level diagrams appear in Figure]^ Each W 
subsequently decays either hadronically (qq), or leptonically {iu, i = e, fj., or r). 
There are then three possible four-fermion final states, hadronic {W~^W~ ^qqqq) , 
semi-leptonic {W~^W~^qqti7), and leptonic {W'^W~^iuiv), with branching 
fractions of 46%, 44%, and 10% respectively. The W^W~ production cross sec- 
tion varies from 3.8 pb at ^/s = 161 GeV to 17.4 pb at = 200 GeV. This can 
be contrasted with the production cross sections for the dominant background 
processes]], a (e+e" Z/7 qq) ~ 100 — 150 pb, a{e'^e~ WeV) ~ 0.6 pb, 
a (e"^e~ Z^e~^e~) ~ 2 — 3 pb, and a (e+e" — > Z/'yZ/'j) ~ 0.5 — 1.5 pb, where 
the spread accounts for variations across the different LEP2 center-of-mass ener- 
gies [45|. The algorithms used to select candidate events exploit the kinematic 
properties unique to the W~^W~ final states. The selection algorithms are sen- 
sitive to all possible W^W~ final states and obtain efficiencies of better than 
about 70% with purities in excess of about 



3.1 Measurement Techniques 

There are two main methods available to measure My/ at LEP2. The first exploits 
the fact that the W^W~ production cross section is particularly sensitive to Mw 
at y/s ~ 2M\\r. In this threshold region, assuming SM couplings and production 
mechanisms, a measure of the production cross section yields a measure of My/. 
In early 1996 each of the LEP experiments collected roughly 10 pb~^ of data at 
y/s = 161 GeV and determined Myy using the threshold technique 

The second method uses the shape of the reconstructed invariant mass distri- 
bution to measure Myy. This method is particularly useful for > 170 GeV 
where the W^W~ production cross section is larger and phase-space effects on the 

®"LEP1" refers to data taken from 1989-1995, when the LEP colhder operated at about 
7s = Afz, while "LEP2" refers to data taken from 1996-2000 at = 161 - 205 GeV. 

^Some of the cross sections given here include kinematic cuts which restrict the final states 
to VF"'"VF~-like parts of phase-space. These cuts are detailed in Reference [la]. 
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reconstructed mass distribution are smaller 
10pb~^ at = 172 GeV |§ in later 1996, 55pb 
180 pb"^ at ^ = 189 GeV in 1998 pi, and 225 pb 



Each experiment collected roughly 

^ at ^/i = 183 GeV in 1997 ~ 

-1 



's 

at ^/s = 192 - 202 GeV in 
1999. Since most of the LEP2 data have been collected at center-of-mass energies 
well above the W~^W~ threshold, the LEP2 M\y determination is dominated by 
the direct reconstruction method. The results reported in this article only include 
the data taken through the end of 1998. 

Each method shall be described in greater detail below. 



3.2 Threshold Determination of Mw 

At center-of-mass energies very near 2 Mw the W~^W~ production cross section, 
aww = c {e~^e~ — > W^W~), is a strong function of M]y, so that a measurement 
of aww can be used to determine the VF-boson mass. This is illustrated in 
Figure which plots the VF-pair production cross section as a function of center- 
of-mass energy for various assumed values of the T^-boson mass. Note that 
for ^/s significantly above or below 2M\y, the various curves converge, so that 
aww has little sensitivity to Mw at those energies. It is only in the threshold 
region that the curves significantly separate, so that a measure of aww affords 
a determination of Mw- To measure the M^-boson mass using the threshold 
method, one needs to 1) select events, 2) determine aww, and 3) extract Mw 
from the aww determination. In practice, steps 2) and 3) are not completely 
independent due to quantum interference effects, which require that special care 
be taken when defining the W~^W~ production cross section. Before discussing 
each of the steps in more detail below, it is useful to detail the aww calculation. 

The V7-pair production cross section is defined to be the production cross 
section for the diagrams given in Figure ^. The separation between the signal 
W~^W~ production diagrams and background diagrams resulting in the same 
four-fermion final states {e.g. e+e" W^W~ — > udud and e"'"e~ — > Z* Z* 
uudd) is complicated by quantum interference effects. In addition, the shape of 
the aww vs. ^/s curve is affected by higher order electroweak and QCD cor- 
rections. The effects of four-fermion interference and the electroweak and QCD 
corrections must all be sufficiently understood theoretically in order to meaning- 
fully determine Mw ■ It turns out that the interference effects are small and can be 
sufficiently addressed in the background subtraction, as discussed in section p. 2. 5 



The effects of the higher-order corrections are larger and will be further discussed 
here. 

As illustrated in Figure |9| [^, the aww vs. ^/s curve gets smeared out near 
2Mw by the effects of the natural width of the VF-boson, Tw, and by the effects 
of initial state radiation (ISR). There are also large corrections associated with 
Coulomb interactions between the two Ws, and some QCD corrections affecting 
W^W~ — > qqf f final states. The theoretical uncertainties associated with cal- 
culating the necessary corrections contribute a 2% uncertainty in aww- These 
theoretical uncertainties are dominated by the uncertainty in the Higgs boson 
mass, which contributes to higher-order electroweak loop corrections and is most 
pronounced near threshold (1.5%). The remaining uncertainties contribute below 



the 0.5% level [SO]. 
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M^=79.83 GeV 
M =80.33 GeV 
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Figure 8: The e+e 
various Mw values. 
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Figure 9: The e~^e~ — > W~^W~ cross section as a function of ^/s using the 
on-shell (Born) approximation, and then including various corrections due to the 
effect of the W^-boson width iTw), the Coulomb interaction between M^+VF", 
and ISRliF 
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3.2.1 EVENT SELECTIONS 

The statistical uncertainty in M^y determined from the threshold method can be 
expressed as 



AMi4/(stat) = ^/<Tww 



dMw 



daww 



(11) 



where eww and V are the W^W" selection efficiency and purity, respectively and 
C is the total integrated luminosity. From this it is obvious that high efficiency, 
high purity selections are important. Separate selections are developed for each 
of the main four-fermion final states — the fully hadronic, the semi-leptonic, and 
the fully leptonic. Each will be discussed separately. The algorithms employed 
are quite involved and vary in the details of their implementation across the LEP 
experiments. In the descriptions below an effort is made to simply emphasize the 
most important discriminating variables and the dominant systematic uncertain- 
ties. For detailed descriptions of the selection algorithms, the reader is referred 
to References [E6|. 



3.2.2 W+W-^qqqq EVENT SELECTION 

The fully hadronic selection is designed to efficiently select W'^W~^qqqq events, 
which are characterized by four (or more) energetic hadronic jets, with little miss- 
ing energy or momentum. The dominant background is from the QCD processes 
e~^e~ — > Z/7 — > qq{+ng), which radiate little energy to ISR. Discrimination re- 
lies primarily on the fact that the jets in signal events tend to be higher energy 
and more spherically distributed than those in background events. In addition, 
to further reduce the QCD background, a kinematic fit can be employed which 
requires the two di-jet masses to be approximately equal. The selections usually 
require high multiplicity, full energy events and exploit the unique W~^W~^qqqq 
kinematics in a multivariate discriminant {e.g. a neural net output) to separate 
signal from background. The typical selection efficiency is about 55% with 80% 
purity. 

For the background estimate, the dominant systematic uncertainty (5%) is as- 
sociated with modeling the dominant QCD background, estimated by comparing 
data to Monte Carlo using high statistics samples of e~^e~ ^ qq events from 

LEPl, and by comparing the estimates from various Monte Carlo generators {i.e. 
Pythia and Herwig). The uncertainty associated with the signal efficiency is 
dominated by comparisons of different generators (Pythia, Herwig, Ariadne, 
KORALW and Excalibur). Contributions from uncertainties in the LEP beam 
energy, E^j^, and ISR are negligible (< 1%). The selection efficiency is also 
negligibly dependent on Mw and on the details of modeling color-reconnection 
and Bose-Einstein correlation effects, which are discussed in more detail in Sec- 
tion Islsl. 



3.2.3 W+W-^qqeu EVENT SELECTION 

The semi-leptonic selection is designed to efficiently select W~^W~ ^qqtu events 
and is typically broken into three separate selections, one for each lepton flavor. 

The W'^W~ ^qqeU and W^W~ ^qq^v events are characterized by two ener- 
getic hadronic jets, a high energy, isolated lepton and large missing momentum 
associated with the prompt neutrino from the leptonically decaying W . The 
dominant background is from radiative e"'"e~ — > Z/^ ^ qq events in which a 
hadron or initial state photon is misidentified as a lepton. Other background 
sources include e'^e~ WeV, e~^e~ ZZ and e~^e~ Ze'^e~ events. The 
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selections require two hadronic jets, an identified, energetic {e.g. E > 25 GeV), 
isolated electron or muon and large missing momentum. The backgrounds from 
radiative e~^e~ qq and e+e" WeV events tend to produce missing momen- 
tum along the beam axis. Requiring a significant missing transverse momentum 
dramatically reduces these backgrounds. The typical selection efficiency is about 
70 — 80% with purities of around 95%. These selections also select about 5% of 
W^W~ —*qqTV events due to the leptonic decays of the tau. 

The dominant systematic uncertainty associated with the selection efficiencies 
(2%) is due to uncertainties in the Monte Carlo modeling of the data and from 
comparing different Monte Carlo generators for e^e~ — > W^W~ events. The 
dominant systematic uncertainty (30 — 50%) associated with the background 
estimate is again due to the modeling of the dominant e^e'^ ^ Z ^ qq back- 
ground and from comparisons of different Monte Carlo generators. For the 
W~^W~ ^qqeV channel, uncertainties from the modeling of four-fermion inter- 
ference, particularly from e+e^ WeV events, estimated by comparing the 
results of different Monte Carlo generators, can increase the total background 
uncertainty to 100%. Since the selections are so pure, these relatively large un- 
certainties in the accepted background cross sections translate into very small 
uncertainties in ctww- 

The W~^W~ -^qqrV events are characterized by two hadronic jets, a r-decay jet, 
and missing momentum associated with two or more neutrinos. The dominant 
background arises from radiative e'^e~ — > Z/7 qq events where a third jet, 
often due to soft gluon emission, is misidentified as a r-jet. These selections are 
very similar to the W^W'^qqeV and W^W~ ^qqiiV selections except that they 
identify the r as a low-mass, low multiplicity (1- or 3-prong), isolated jet. Since 
the lepton identification is looser than for the other W~^W~ —>qq£u selections, the 
background tends to be higher. Selection efficiencies vary widely among the LEP 
experiments — from 35 — 45% exclusive of those W~^W~ ^qqrV events identified 
by one of the other selections. Due to the looser lepton identification requirements 
these algorithms typically select an additional 3 — 5% of W^W~ ^qqeV and 
W'^W~ ^qqjiV events failing the above selections. The typical purity of this 
selection also varies widely among the LEP experiments, 65 — 85%. 

The dominant systematic uncertainty associated with the selection efficiency 
(2.5%) is due to the modeling of the lepton identification variables, estimated 
by comparing LEPl data and Monte Carlo, and from the comparison of vari- 
ous Monte Carlo generators. The dominant systematic uncertainty associated 
with the estimate of the accepted background cross section (20%) comes from 
the modeling of hadron mis-identification, estimated by comparing the data and 
Monte Carlo fake rates in LEPl e^e^ —f qq events. 

3.2.4 W+W~^eueu EVENT SELECTION 

The fully leptonic channel, l^+VF"— >£f£f, is characterized by two high energy, 
isolated, acoplanar leptons. The selections typically start by requiring a low 
multiplicity and large missing transverse momentum. There are six distinct 
final states (ee, efj,, er, /x//, /xr and rr), which have differing dominant back- 
ground sources. Potential background sources are two photon, e+e" WeV, 
e'^e~ — > Z^e'^e~ , and radiative e'^e~ Z^ ^ £^£~ events. In general the 
W^W'^tPtP selection involves several independent and overlapping sets of cuts 
which employ various combinations of specific electron, muon, and tau identifi- 
cation algorithms. Backgrounds are usually rejected by requiring large missing 
energy, large transverse momentum and a large lepton-lepton acoplanarity. The 
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efficiency varies widely across the LEP experiments, from about 45% for Delphi 
and L3, to about 65% for Opal and Aleph. The selection purities are around 
90 - 95%. 

The dominant systematic uncertainty associated with estimating the selection 
efficiency (2%) is due to the modeling of lepton identification variables, specif- 
ically those sensitive to FSR modeling {e.g. isolation variables). Comparisons 
from different Monte Carlo generators also contribute. The dominant systematic 
uncertainty associated with estimating the accepted background cross sections 
(70%) arises from limited Monte Carlo statistics and comparisons of different 
Monte Carlo generators. The effect of detector mis-modeling is small owing to 
the experience gained at LEPl using e~^e~ ^ £~^£~ data. However, veto 

cuts employed in these low multiplicity final states are particularly sensitive to 
beam related backgrounds which are not included in the Monte Carlo. These are 
estimated using random trigger crossings and have the consequence of reducing 
both the signal efficiency and accepted background by a factor of order 0.95 — 1.0 
with a relative uncertainty of a few percent. 



3.2.5 DETERMINING aww AND Mw 

A maximum likelihood procedure is used to determine ctww- The likelihood is 
usually taken to be the product of Poisson probabilities of observing Ni events 
when expecting fii{aww) = ^■WwiV^i-^i+o'i,'^^] events, where £i, Bi, and crgg^. 
are the selection efficiency, branching ratio, and accepted background cross sec- 
tion, respectively, for the ith. selection, and C is the integrated luminosity. In the 
likelihood calculation, correlations between the channels are properly accounted 
for and Standard Model branching ratios are assumed. The accepted background 
cross section is assumed independent of M-^. Four-fermion interference effects 
are typically accounted for either by applying a correction factor, fi, to the prod- 
uct aww^i (^ww^ifi, or by adding a factor, Fi, to a?ff^ a?'i'^, + Fj. These 
corrections are negligible (|1 — /j| < 1%) for all channels except the qqeU and 
iViU channels which have fqqev ^ ^-^^ i^qqeV ~ —0.05 pb) and ~ 0.97 

(-^iuiu ~ +0.01 pb). These correction factors are determined by comparing the 
predicted total accepted cross section {i.e. signal plus background) calculated 
including and excluding these interference effects. A systematic uncertainty is 
estimated by comparing the predictions from different Monte Carlo generators 
which include the four-fermion (4f) interference effects. Given the large statisti- 
cal uncertainty on aww these corrections do not significantly affect the result. 
Strictly speaking, the correct manner in which to account for these interference 
effects, is to make the replacement aww ~^ o'4f in the likelihood calculation, 
where a^f is the M^y-dependent 4-fermion production cross section, including 
the interference effects. 

The measured W~^W^ production cross section is then compared to SM predic- 
tions for aww dependent on Mw and Ebm- The likelihood equation is modified 
so that aww '^ww{Mw, Ebm). The dependence of aww on Mw is calculated 



using the semi-analytic program Gentle |51], which includes the higher order 
electroweak and QCD corrections. The results from each of the LEP experiments 
are given m Table |. 

In addition to the systematic uncertainties associated with the selection ef- 
ficiencies and accepted background cross sections described above, there are 
uncertainties due to higher-order corrections affecting aww{Mw,Ebm) (2% at 
^/s = 161 GeV as discussed above) and to the precision of the LEP determina- 
tion of Ehm (±27 MeV at = 161 GeV |1). 
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Table 3: I^-pair production cross section, aww, and H^-boson mass, Mw, 
results for data taken at ^/s = 161 GeV 



exp 


aww ± (sta) lb (sys) (pb) 


Mw ± (sta) ± (sys) (GeV) 


£ (pb-^) 


Aleph 


4.23 ±0.73 ±0.19 


80.14 ±0.34 ±0.09 


11 


Delphi 


3.67l:[!|^ ±0.19 
2.89l[J^o ±0-14 
3.62t|]|^ ±0.16 


80.40 ± 0.44 ± 0.09 


10 


L3 
Opal 


sn on+0.47+0.09 

OU.»U_o.4i_o.08 

80.40l|];^f ± 0.10 


11 
10 



3.2.6 COMBINATION OF Mw DETERMINATIONS FROM THRESHOLD 

The combined determination of Mw is extracted from the LEP combined mea- 
surement of aww, which is then compared to the Gentle prediction for awwi^w , Etm) 
assuming SM couplings and the LEP average center-of-mass energy as shown in 
Figure When making the combination the expected statistical uncertainty is 
used to calculate the weight of each experimental measurement so as to avoid 
biasing the result. Since the statistical uncertainties dominate each individual 
measurement as well as the combined result, the smallest quoted systematic un- 
certainty (0.14 pb) is conservatively taken to be fully correlated between exper- 
iments. Note that since an individual experiment's weight in the combination 
is driven by its statistical uncertainty, this procedure does not affect the cen- 
tral value of the combination and yields a conservative estimate of the combined 
systematic uncertainty. Combining the four LEP experiments' determinations of 
o'ww yields 

aww = 3.69 ± 0.45 pb (12) 

with a pe^ degree of freedom of 1.3/3. Using the LEP average center-of-mass 
energy of 161.33 ± 0.05 GeV and the Gentle prediction, the VK-boson mass is 
then determined from this threshold (TH) method to be | |53| 

Mw'(TH) = 80.400 ± 0.220(exp) ± 0.025{Ebm) GeV. (13) 

The statistical uncertainty dominates the experimental uncertainty, which has 
a contribution of approximately 70 MeV from correlated systematics. If this 
method were to be employed in the future {e.g. at an NLC), a potentially lim- 
iting uncertainty is due to the modeling of fragmentation and hadronization, 
which has a large effect on the W^W~ ^qqqq channel and is correlated among 
the experiments. This uncertainty is presently the single largest contribution to 
the total uncertainty assigned to the LEP combined aww at higher center-of- 
mass energies (where the statistical uncertainties are smaller) and contributes 
approximately a 50 MeV uncertainty to Mvt'(TH). 



3.3 Direct Reconstruction of Mw 

As demonstrated in Figure |^, at center-of-mass energies above 170 GeV the 
W^W~ production cross section becomes significantly less sensitive to Mw- At 
these energies one can extract a measure of Mw from the invariant mass dis- 
tribution of the W^-decay products. The sensitivity to uncertainties associated 
with the modeling of events near the phase-space limit {Mw+ = Mw- = E hm) is 
greatly reduced since {^/s — 2Mw) ^ ^w- However, as discussed in Section 3.3.4 , 
the modeling of various final state interactions becomes more important in the 
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Figure 10: The LEP combined aww (shaded band) is compared to a semi- 
analytic calculation of awwi^w, \/^) using the LEP average center-of-mass 
energy to extract a LEP combined determination of My\r using the threshold 
method [El. 



W~^W~^qqqq channel. To measure Mw using this direct reconstruction tech- 
nique one must 1) select events, 2) obtain the reconstructed invariant mass, rrirec, 
for each event, and 3) extract a measure of Mw from the rrirec distribution. Each 
of these steps is discussed in detail in the sections below. 



3.3.1 EVENT SELECTION 

The expected statistical uncertainty on Mw determined from direct reconstruc- 
tion will vary inversely with the selection efficiency and purity. At higher center- 
of-mass energies the l^+Vl^" production cross section increases by over a factor 
of four, while the dominant background cross sections increase less rapidly, or 
even decrease. This affords selections with greater efficiencies for the same puri- 
ties relative to those employed for the selection at ^/s = 161 GeV. Nevertheless, 
the algorithms employed are very similar to those described in Section ^.2.1| and 
so will not be further discussed here. Typical selection efficiencies (purities) are 
85% (80%) for the W+W-^qqqq channel, 90% (95%) for the W+W-^qqeV and 
W~^W~ ^qqfiu channels, and 65% (85%) for the W~^W~ ^qqrV channel. The 
W~^W~^ii'£T' channel does not significantly contribute to the determination of 
Mw from direct reconstruction and is not discussed further. For the high energy 
data taken through 1999, these efficiencies and purities yield approximately 7000 
W~^W~ — > qqf f candidate events, about 1100 of which are Tioii-W~^W~ back- 
ground. The selection efficiencies have a total uncertainty of about 1% (absolute) 
and have a negligible effect (< 1 MeV) on the determination. The accepted 
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background cross sections have a total uncertainty of 10 — 20% (relative) and 
affect the M\y determination at the 10 — 15 MeV level (cf. Section |3.3.4D . 

3.3.2 INVARIANT MASS RECONSTRUCTION 

For each selected event, an invariant mass is reconstructed from the W decay 
products. There are several methods available for reconstructing the invariant 
mass oi aW candidate. The best resolution is obtained by using a kinematic fit 
which exploits the fact that the center-of-mass energy of the collision is known a 
priori^ Since the type of fit employed varies for each final state, each will be dis- 
cussed separately. While the details of the fits differ among the LEP experiments, 
the important features are similar. 

Selected W~^W~ ^qqqq events are forced into a four-jet configuration using, 
for example, the Durham algorithm |54]. A kinematic fit is then performed 



to estimate the reconstructed invariant mass of the W candidates. A fit which 
incorporates the constraints of energy and momentum conservation (4C fit) yields 
two reconstructed invariant masses per event {mreci,iTT'rec2) i one for each 1^-boson 
in the final state. A fifth constraint can be incorporated by neglecting the finite 
W width and constraining the invariant masses to be equal, nireci = '^rec2- For 
each event, this 5C fit yields a single reconstructed mass, rrirec: its uncertainty, 
CTree ! and a fit x^-probability. The fit requires as input the jet momenta, energy, 
and their associated uncertainties. A complication of the qqqq final state is due 
to the pairing ambiguity — there exist three possible jet-jet pairings for a four- 
jet final state. This pairing ambiguity gives rise to a combinatoric background 
unique to this channel. Each LEP experiment employs a different technique for 
differentiating among the combinations. Typically an experiment will use the 
best one or two combinations. The correct combination is among those used in 
about 85 — 90% of the events. For events with the correct pairing, the kinematic 
fit has a resolution of about 0.7 GeV per event and is dominated by the jet angular 
resolution. The wrong combinations are treated as a background. It should be 



noted that the shape of the combinatoric background is fairly flat (cf. Figure 11). 
Because of this the M\y determination is not critically dependent on how well 
known the fraction of correct pairings is. 

Selected qqeU and qqfJl^ events are forced, after removing the lepton candi- 
date, into a two jet configuration. All four LEP experiments use a kinematic fit 
employing energy and momentum conservation constraints and the equal mass 
constraint. Since the prompt neutrino from the leptonically decaying W takes 
three degrees of freedom, this is a 2C fit yielding a single reconstructed mass, 
uncertainty and fit x^'Pi'obability per event. The fit requires as input the jet 
and lepton energy and momenta and their associated uncertainty. The qqel7 and 
qqfu' events have a resolution of roughly 1.0 GeV and 1.1 GeV, respectively, per 
event. This resolution is dominated by the uncertainty in the lepton energy. 

Selected qqrV events are forced, after removing tracks and clusters associated 
with the r-decay, into a two jet configuration. The treatment of qqTl7 events 
varies among the LEP experiments. Basically the invariant mass of the hadronic 
system is used, the resolution of which can be improved by requiring energy and 
momentum conservation and employing the equal mass constraint. The resolution 
of the qqrV events is approximately 1.5 GeV per event and is dominated by the 
resolution of the jet energies. 

^Strictly speaking, this is not true since any ISR reduces the coUision energy to less than 
twice the beam energy. The kinematic fits assume no ISR. The effe ct of ISR uncertainties is 



incorporated in the total systematic uncertainty discussed in Section 3.3.4 and is small. 
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3.3.3 EXTRACTING Mw 

The ensemble of selected events yields a rrirec distribution from which a measure of 
Mw is extracted. There are several methods available for extracting Mw- Aleph, 
L3, and Opal all employ a traditional maximum likelihood comparison of data to 
Monte Carlo spectra corresponding to various Mw- In addition to its simplicity, 
this method has the advantage that all biases {e.g- from resolution, ISR, selection, 
etc.) are implicitly included in the Monte Carlo spectra. The disadvantage of this 
method is that it may not make optimal use of all available information. Delphi 
employs a convolution technique, which makes use of all available information; 
in particular, events with large fit-errors are de-weighted relative to fits with 
small fit-errors. The convolution has the limitations that it requires various 
approximations (e.g. the resolution is often assumed to be Gaussian) and often 
requires an a posteriori correction as the fit procedure does not account for all 
biases, notably from ISR and selection. All experiments employ an analytic fit 
of a relativistic Breit-Wigner (with s-dependent width) -l-background to the data, 
which also requires a posteriori corrections, as a cross-check. Since their dominant 
systematic uncertainties differ, Mw is measured from the qqqq and the qqiP 
samples separately. These are then combined, taking into account correlations, 
to yield a final measurement of Mw- In the results given here, the Standard 
Model relation between Mw and Tw has been assumed |5C]. 



Table 4: LEP Results for the qqtP channel for data taken at ^/s = 172 — 189 GeV 



exp 


Mw ± (stat) ± (syst)/GeV 


Aleph 


80.343 ± 0.089 ± 0.040 


Delphi 


80.297 ±0.141 ±0.064 


L3 


80.224 ±0.117 ±0.067 


Opal 


80.362 ± 0.090 ± 0.053 


LEP 


80.313 ±0.052 ±0.036 



Table 5: LEP Results for the qqqq channel for data taken at ^/s = 172 — 189 GeV 



exp 


Mw ± (stat) ± (syst) ± (CR/BE)/GeV 


Aleph 


80.561 ± 0.095 ± 0.050 ± 0.056 


Delphi 


80.367 ± 0.094 ± 0.037 ± 0.054 


L3 


80.656 ± 0.104 ± 0.071 ± 0.092 


Opal 


80.345 ± 0.098 ± 0.074 ± 0.055 


LEP 


80.429 ± 0.049 ± 0.046 ± 0.058 



The results from each LEP experiment, using data collected at y/s = 172 — 
189 GeV 0, m, H], are given in Table | for the qq&7 channel and in Table | 
for the qqqq channe^. Also included is the mass obtained when combining all 

^These results are based in part on preliminary numbers for the data taken at -^s — 189 GeV. 
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the measurements. In the combination correlations are taken into account as 
described in Section 3.3.5 . Figure shows the Opal fits for the data taken at 
^/s = 189 GeV. 

3.3.4 SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties for a "typical" LEP experiment are given in Ta- 
ble ^. This table should be taken as a general guide. The actual numbers vary 
by as much as ±20 MeV from experiment to experiment. It is still the case 
that the total uncertainty of a single experiment is dominated by the statistical 
uncertainty. As this is the case, the experiments are at various stages in develop- 
ing more sophisticated methods to estimate the limiting systematic uncertainties. 
This largely accounts for the range of uncertainties and not any inherent detector 
or methodological advantages. It should be noted that for all four LEP exper- 
iments the uncertainties associated with ISR, fragmentation, and four-fermion 
interference are limited by the statistics of the comparison. Uncertainties associ- 
ated with the selection efficiencies and accepted backgrounds are included in the 
line labeled "fit procedure" . For the qqiP channel the largest single contribution 
to the systematic uncertainty is due to detector effects {e.g. energy scales, res- 
olutions, and modeling). These uncertainties are expected to decrease as more 
data are collected. For the qqqq channel the dominant systematic uncertainty is 
due to color-reconnection and Bose-Einstein correlation effects. Each source of 
uncertainty and the methods for estimating it are briefly described below. 

Initial state radiation: Uncertainties due to the modeling of ISR are estimated 
by comparing the m^ec distributions of Monte Carlo samples which include 
ISR corrections to differing orders. The standard Monte Carlo include cor- 
rections to next-to- leading- log O(a^). The differences are negligible in sam- 
ples of several million events. The uncertainty is conservatively taken to be 
10 MeV. 

Four-fermion interference: The systematic uncertainty associated with the 
modeling of four-fermion interference effects is usually estimated by com- 
paring fit results using Monte Carlo reweighting techniques which employ 
matrix element calculations including and excluding the interferences. The 
differences are small and the comparisons often statistically limited. 

Fragmentation: A variety of methods is employed to estimate the uncertainty 
associated with fragmentation modeling. Typically, LEPl data are used to 
constrain model parameters. Those parameters are then varied in several 
Monte Carlo samples, which are re-fit as data. The results are compared 
to a Monte Carlo sample with default parameters. The differences are typ- 
ically small except for the W~^W~^qqqq channel, where they are on the 
order of 30 MeV. Since this systematic is strongly correlated among the 
LEP experiments, it is one of the dominant systematic uncertainties in the 
LEP combined measurement. Delphi employs an alternative method for 
estimating these uncertainties |48], in which two — > qq events are boosted 



to the appropriate center-of-mass energy and overlayed in data and Monte 
Carlo; the two samples are treated as data and their resulting fit masses 
compared. 

Detector effects: Calibration data collected at the start and end of each year at 
^/s = Mz, are used to establish limits on the uncertainties associated with 
detector modeling in the Monte Carlo. Uncertainties in energy scale and 
resolution are estimated using — > e'^ e~ , fi'^ ^~ , qq data. The linearity is 
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Figure 11: The preliminary L3 fit result for the W^W 
data taken at y/s = 189 GeV. 
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Figure 12: Preliminary Opal fit results using data collected at y/s = 189 GeV. 
The points are the data and the open histogram is the fit result. The non-WW 
background contribution, as estimated from Monte Carlo, is shown as the hatched 
histogram. 
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checked using Bhabha, e+e~7, and 3-jet events in data collected at higher 
center-of-mass energies. The angular resolutions are similarly determined. 
Typically the jet and lepton energy scales are known to within 0.5%, while 
the modeling of the angular and energy resolutions have associated uncer- 
tainties of the order 5 — 20% depending on polar angle. These uncertainties 
are propagated to the Mw measurement by comparing the fit results of 
Monte Carlo samples in which the appropriate quantity has been scaled or 
varied to the results from the default Monte Carlo sample. The observed 
differences are used as estimates of the associated uncertainty on M]y due 
to detector modeling effects. 

Fit procedure: The type of effects considered depend on the fit method em- 
ployed to extract Mw from the Tfi^QQ distribution. These include uncertain- 
ties associated with the background normalization and shape, and fit biases. 
In general the total effect is very similar across methods. The background 
normalization is varied within uncertainties determined from a dedicated 
(Jww analysis yielding small effects. The shape of the background nirec 
distribution is cross-checked using data where possible and compared across 
different Monte Carlo generators otherwise, also yielding small effects. The 
linearity of the fit methods is determined from Monte Carlo samples gen- 
erated assuming various My/ values. These samples are also used to verify 
that the statistical uncertainty is accurate. For the convolution and Breit- 
Wigner methods, these Monte Carlo samples are used to calculate the nec- 
essary bias corrections, whose uncertainties are then propagated to the final 
uncertainty. 

Beam energy: The uncertainty in the beam energy enters via the constraints 
employed by the kinematic fit and should be of order dMw = Mw . The 
effect on the measured M\y is estimated by re-fitting all data changing the 
£'bm used in the fit and calculating the mean difference in fitted rrirec on an 
event-by-event basis. The beam energy is estimated using the method of res- 
onant depolarization |^5|; the last energy point for which the depolarization 
method works is Ei^^ ~ 60 GeV. An extrapolation is required to estimate 
the beam energies at which the data are taken, Ei,^ « 100 GeV. The un- 
certainty in the beam energy is about 20 MeV and is dominated by the 
uncertainty in this extrapolation. With the addition of more resonant de- 
polarization data and new techniques, it is expected that the uncertainty on 



M^y due to uncertainties in E'^.^ will be reduced to around 10 MeV |52, 
The spread in the beam energy, cr^^ ~ 150 MeV [^], has a negligible 
effect on My/. 

CR/BE: Since the typical decay distance of the VF-bosons, l/F^y 0.1 fm, is 
much smaller than the typical fragmentation radius, 1/Xqcd ~ 1 f™, the 
decay products originating from different VF-bosons cannot be considered 
as independent — i. e. they can "talk" to each other. The modeling of this 
cross-talk in the Monte Carlo spectra used to extract Mw is an additional 
source of systematic uncertainty in the W^W~ ^qqqq channel. The cross- 
talk can arise through two mechanisms, Bose-Einstein correlations (BE) 
and color reconnection effects (CR) [^][|^. The modeling uncertainty is 
estimated separately for BE and CR and is model dependent in both cases. 
In each case Monte Carlo samples employing implementations of various 
CR/BE models are treated as data and an Myr bias is estimated. The 
systematic uncertainty is chosen to include the full range of variation among 
the models explored and is on the order of 50 — 60 MeV. 



30 



Glenzinski & Heintz 



It should be noted that there has been recent progress in experimentally con- 
straining the available CR models by comparing event shape and charged parti- 
cle multiplicity distributions as predicted by various Monte Carlo models (both 
including and excluding CR effects) with those observed in the high-energy data. 
In addition, studies using LEPl data can also be used to test the available mod- 
els 1 58]. On the basis of these studies, some of the models have been excluded 
as they fail to adequately describe the data [59|, thus enabling a reduction in 
the associated systematic uncertainty (from ~ 100 to ~ 50 MeV). For a more 
complete discussion, see Reference ||5^. Additional data should help to further 
constrain the remaining CR models and thus improve this uncertainty further. 



Table 6: Table of systematic uncertainties on Myy from direct reconstruction 
for a "typical" LEP experiment 



systematic 


uncertainty (MeV) 


source 


qqtP 


qqqq 


initial state radiation 


10 


10 


four-fermion 


10 


10 


fragmentation 


25 


30 


detector effects 


30 


30 


fit procedure 


20 


20 


Sub-total 


46 


49 


beam energy 


17 


17 


CR/BE 




60 


Total 


49 


79 



3.3.5 COMBINATION OF Mw DETERMINATIONS FROM DIRECT RE- 
CONSTRUCTION 

Each of the LEP experiments provides their measured VF-boson mass for the fully- 
hadronic and semi-leptonic channels separately for each center-of-mass energy 
along with a matrix of associated uncertainties. The uncertainties are broken 
down into four components: 

1. uncertainties uncorrelated between channels and experiments {e.g. the sta- 
tistical uncertainty or background normalization and shape uncertainties) 

2. uncertainties correlated among the channels of a given experiment, but un- 
correlated between experiments {e.g. detector modeling uncertainties) 

3. uncertainties uncorrelated between the channels, but correlated among the 
experiments {i.e. CR/BE uncertainties) 

4. uncertainties correlated between the channels and among the experiments 
{e.g. ISR, fragmentation, Ebm uncertainties). 

In this way the correlations between channels and among experiments are ac- 
counted for. The correlation of the Ebm uncertainty across the different years 
is also taken into account. The results for the combined qqqq and qqlV chan- 
nels are given in the last lines of Tables |^ and || and are 25% correlated with a 
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X^/dof = 17.9/20. Combining all the direct reconstruction (DR) results into a 
single mass yields 

Mw (DR) = 80.347 ± 0.036(stat ) ± 0.036(syst ) ± 0.020(CR/BE) ± 0.017(^b„) GeV 

where the uncertainties associated with CR/BE modeling and with the LEP beam 
energy have been listed separately ||6^. The dominant systematic uncertainty is 
associated with the fragmentation model, which is correlated among the experi- 
ments (they all employ the same models in their Monte Carlo) and contributes 
an uncertainty of approximately 20 MeV. The effect of the CR/BE uncertainty 
is to de-weight the qqqq measurements relative to the measurements in the qqiu 
channels. 



3.4 Combination of LEP Results 

The Mw determination using the threshold method is combined with the determi- 
nation using the direct reconstruction method taking account of the correlations. 
In particular, the systematic uncertainties associated with the LEP beam energy, 
and the modeling of ISR, fragmentation, and four-fermion interference effects are 
taken as correlated. Note that the weight of the threshold determination of M\y 
in the combination is driven by the statistical uncertainty of that measurement. 
The LEP combined result, assuming the Standard Model relation between the 
W decay width and mass, is 

Mvi/(LEP) = 80.350 ± 0.056 GeV, (14) 

where the uncertainty includes both statistical and systematic uncertainties and 



is dominated by the determinations using direct reconstruction methods [ 30 



4 WHAT DO THESE MEASUREMENTS TELL US? 



4.1 Combination of Results 

Direct measurements of the VF-boson mass have been performed in two kinds of 
experiments, the production of T^-bosons in pp collisions and the production of 
W~^W~ pairs in e~^e~ collisions. 

Until 1996, pp collisions were the only source of PF-bosons. The advantage of 
pp colliders lies in the large W production cross section and the low background 
levels. The pp data give about 100,000 W ^ tu candidate events with about 97% 
purity. The production of Z-bosons, dynamically and kinematically very similar 
to VF-boson production, provides a very convenient control data sample. The 
disadvantage of pp colliders is that the parton center-of-mass frame is not known 
on an event-by-event basis and therefore systematic effects arising from the struc- 
ture of the protons must be understood. The combined pp-collider measurement 
is Mwipp) = 80.452 ± 0.060 GeV. 

Since 1996, e~^e~ collisions with enough energy to produce pairs of VF-bosons 
are available. The advantage of e~^e~ collisions is that since the initial particles 
are point like and so the center-of-mass energy of the collision is known, kinematic 
fits can be employed to fully reconstruct events and thus yield invariant mass 
resolutions comparable to the T^-boson width. The disadvantage of e~^e~ colliders 
is that the T^"'"M^~ production cross section is two orders of magnitude smaller 
than at pp colliders, resulting in smaller and less pure event samples (about 
22,000 events with about 90% purity). In addition, the modeling of final-state 
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interactions in W~^W -^qqqq events must be understood. The combined LEP 
measurement is Miy(LEP) = 80.350 ± 0.056 GeV. 

The two determinations of the VF-boson mass are completely uncorrelated. A 
combination of both results is simple, resulting in a world average of 

Mw = 80.398 ± 0.041 GeV (15) 

with a of 1.6. Having two independent, precise determinations of this param- 
eter in agreement with each other lends significant credibility to the results. 

Within the framework of the Standard Model, the measurement of the VF-boson 
mass determines the radiative corrections, Ar, in Equation]^. These corrections 
have a large contribution from the running of the electromagnetic coupling. We 
can absorb this into the value of a by writing 



1 - Ar 1 - Ar, 



ewk 



(16) 



For the residual contribution from electroweak loops, we find 
Arg^fe = —0.0268 ± 0.0027, about 10 standard deviations from zero. 

4-2 Comparisons and Constraints within the Standard Model 

The Standard Model provides us with a framework that allows us to relate the 
measurements from many processes that involve the electroweak interaction. The 
main sources of such information are measurements of the properties of the Z- 
boson at LEPl and SLC, the study of deep inelastic neutrino scattering at Fer- 
milab, and the measurement of the mass of the top quark at the Tevatron. 
LEPl and SLC have provided us with a wealth of very precise measurements 



of the properties of the Z-boson [31|. At tree level, the properties of the Z- 
boson are determined by its mass, the weak mixing angle, and the fine structure 
constant. Radiative corrections are dominated by the masses of the top quark 
and the Higgs boson. The wonderful success of the Standard Model lies in all 
measurements being consistent with single values of these parameters. The mass 
of the Z-boson is measured directly from the line shape, and the fine structure 
constant, evolved to = Mz^ is derived from measurements of R, the ratio of 
the e^e~ cross sections to hadrons and to ^'^ ^~ . The other three parameters 
are extracted from a fit to the measurements. The W mass then follows from 
Equation 

The CCFR [|6^ and NuTeV experiments at Fermilab measure the ratio of 
charged current and neutral current interactions of neutrinos. This ratio depends 
directly on 1 — Mw'^/Mz^- From the measured value 

1 - Mw^/Mz^ = 0.2255 ± 0.002l|^, a value for the W"-boson mass of Mw = 
80.250 ± 0.109 GeV can be derived. 

At the loop level many other parameters contribute (mostly negligible) cor- 
rections to the tree level values. Due to the large mass difference between the 
top and bottom quarks, radiative corrections involving top-quark loops are im- 
portant. The CDF and D0 Collaborations have measured the top quark mass 



directly |64]. Their combined value is rritop = 174.3 it 5.1 GeV. 

A fit of the Standard Model to all measurements except the direct measure- 
ments of the VF-boson mass returns [^] M\y = 80.381±0.026 GeV as its preferred 
value. This value is in excellent agreement with the combined direct measure- 
ments, in support of the validity of the Standard Model. The mass of the Higgs 

i"CCFR and NuTeV combined. 
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boson is the only parameter which is not measured experimentahy. Loops con- 
taining Higgs bosons also contribute important radiative corrections. A fit to 
all electroweak data, including the measurements of the VF-boson mass, prefers 
Mh = 77tll GeV for the mass of the Higgs boson [65| ]. 

The Higgs-boson mass can also be constrained based on the measured values of 
the W^-boson and top-quark masses alone. This is shown graphically in Figure ^ 
The shaded bands indicate the values of the T^-boson mass predicted by the 
Standard Model as a function of the top quark mass, for given values of the 
Higgs-boson mass [|6q|. The width of the bands indicates the variation due to the 
uncertainty in a{Mz ) |^^, which dominates the uncertainty in the predictions. 
The ellipse indicates the two-dimensional 68% confidence-level interval defined by 
the measured mass values. The inset shows a plot of between the measured 
values and the predictions as a function of the Higgs-boson mass. The preferred 
Higgs-boson mass is ritfi GeV. Values above 277 GeV are excluded at 90% 
confidence level. The dashed contour shows the 68% confidence-level interval 
from the fit to all other electroweak data [S5|. 




Figure 13: Comparison of direct measurements of the VF-boson and top- 
quark masses with indirect measurements |65] and predictions of the Standard 
Model |6f:]. The indirect constraint is in part based on preliminary data. 
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4.3 Constraints Outside the Standard Model 

Any particle that couples to the VF-boson can contribute loop corrections to the 
value of the VF-boson mass. Thus a measurement of the T^-boson mass does not 
only test the Standard Model, but it is, at least in principle, also sensitive to 
non-standard physics. In the minimal supersymmetric model, corrections could 
increase the W^-boson mass by as much as 250 MeV [^8|. The correction is larger, 
the lower the scale of supersymmetry breaking. If the scale of supersymmetry 
breaking is more than a few hundred GeV, supersymmetry decouples from Stan- 
dard Model physics and the effects of supersymmetric loop corrections on the 
W^-boson mass become small. Supersymmetric particles that give large correc- 
tions must be relatively light and would also be the first ones to be seen in direct 
searches. Thus, precision measurements are unlikely to increase the sensitivity of 
direct searches for supersymmetric extensions of the standard model. 



5 FUTURE PROSPECTS 

The precision on the world average VK-mass is expected to improve significantly 
over the next 5 years, and could possibly improve dramatically over the next 
decade or so. 

By the end of 2000, the LEP experiments will have more than doubled the 
statistics of their W~^W~ data sets relative to what has been included in this re- 
view. The inclusion of the additional data will yield a statistical uncertainty for 
the combined LEP measurement of M\y of about 25 MeV. Already a significant 
effort has been made to reduce the systematic uncertainties, particularly those as- 
sociated with the detector energy scales and resolutions. It is expected that these 
errors can be brought to the 20 MeV level. The additional constraints afforded by 



the LEP spectrometer project |5C] and by additional depolarization data make it 
likely that the uncertainty associated with the LEP beam energy will be reduced 
to roughly 10 MeV. It is difficult to predict how the systematic uncertainties 
associated with the modeling of Bose-Einstein and color-reconnection effects in 
the W~^W~ ^qqqq channel will evolve. While it is true that the additional data 
will provide more stringent tests, it is unknown whether the additional sensitiv- 
ity will actually reveal a discrepancy large enough to reject any of the remaining 
viable models. Assuming none of the remaining models are rejected, so that the 
CR/BE uncertainty remains the same, the projected total uncertainty on Mw at 
the end of LEP2 would be of order 35 — 40 MeV from the LEP combination. If 
the CR/BE uncertainty can be reduced to less than 15 MeV, it may be possible 
for the LEP combined determination to reach a total uncertainty of 30 MeV. 

In spring of 2001, the CDF and D0 experiments will start taking data at 
the Fermilab TeVatron. They anticipate collecting 2 — 3 fb~^ of data by 2004, 
which should give a five-fold improvement in the statistical uncertainty of the 
TeVatron Mw measurement. The systematic uncertainties associated with the 
energy scale and other detector effects are dominated by the statistics of the 
Z control samples and are expected to scale accordingly. On the other hand, 
the systematic uncertainty associated with the H^-production modeling does not 
scale directly with statistics and may improve only moderately to about 20 MeV. 
The uncertainty from the combined TeVatron Mw determination is expected to 



be about 30 MeV||]. 



On the time-scale of the next five years, it is expected that the world aver- 
age M^-mass will have a total uncertainty of 20 — 25 MeV — a factor of two 
improvement relative to the present uncertainty. The SM constraint on Mh af- 
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forded by the Mw measurements alone will be comparable to that afforded by the 
sm9w measurements of LEP and SLD, which presently yield an uncertainty of 
order AMh = Mh . Although these two sets of constraints are correlated through 
^^top-dependent corrections, it is still interesting to compare them since they have 
differing sensitivities to the various radiative correction terms. A significant im- 
provement to the SM constraints on Mh requires a more precise determination 
of rritop- On the same time-scale, the TeVatron experiments are expected to mea- 
sure mtop to 2 — 3 GeV[69]. Including this improvement should yield constraints 
on Mh with uncertainties on the order of AMh = 0.5Mh (assuming that the fits 
continue to give a central value of Mh ~ 0(100) GeV) [|6^] . 

Looking further ahead, on the time-scale of 5-10 years, it is possible that the 
LHC experiments, CMS and ATLAS, will measure the W-mass to a precision 
of 15 MeV and mtop to 2 GeV|7C]. And on the time-scale of > 10 years, future 
high luminosity e^e~ or //^jU~ colliders might yield the statistics to envision a 
< 10 MeV measurement of Myy using the threshold method and a < 1 GeV 
measurement of rritop [ 71 1| 72]. If achieved, such precision measurements would 
yield constraints on Mh with uncertainties of 0{1 — 10) GeV - which we can 
hope will offer, by that time, a very interesting comparison with the directly 
measured M/^! 



6 CONCLUSIONS 

The mass of the W boson has been measured by many experiments at pp and 
e~^e~ colliders. All measurements are in good agreement. The world average of 
all measurements of the VF-boson mass is 80.398 it 0.041 GeV. Based on measure- 
ments of other parameters, the Standard Model of the electroweak interactions 
leads to a prediction of 80.381 ± 0.026 GeV for the mass of the W boson, in 
excellent agreement with the measured value. In the framework of the Standard 
Model this measurement of the VF-boson mass, together with the measurement of 
the top-quark mass, constrains the Higgs-boson mass to values below 280 GeV at 
90% confidence level. Over the coming decade, a reduction in the uncertainty of 
the direct measurement of the VF-boson mass of at least a factor two is expected. 
As the top quark mass is measured more precisely and the reach of searches for 
the Higgs boson increases, the comparison of the indirect constraint on the Higgs- 
boson mass and its direct measurement or exclusion region will become one of 
the most interesting tests of the Standard Model. This test will for the first time 
close in on the symmetry breaking sector of the Standard Model about which 
very little is presently known. 
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